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Abstract. The impact of a typical Shoemaker- Levy 9 fragment produced three 
light peaks as seen from Earth. The first peak is related to the entry of the fragment 
into the Jovian atmosphere. The second peak occurs when the exploding fireball 
rises above Jupiter’s limb into direct view from Earth. The third peak, much the 
brightest, occurs when the ejecta plume falls back on the atmosphere. By contrast, 
Galileo, which had a direct view of the impacts, saw two peaks, one at entry, and 
one at plumefall. Here we present a simple, highly idealized model of a ballistic 
plume, which we then use to fit the observed light curve of the R impact as recorded 
at Mauna Kea and Mount Palomar. From the light curve we find that the nominal 
R fragment had diameter 450-500 m and mass ~2 - 3 X 10 13 g. The uncertainty in 
the mass is probably about a factor of 3, with a smaller event more likely than a 
larger one. 


The Light Curve 

The general light curve for an impact of a fragment of 
P/Shoemaker-Levy 9 (SL9), as viewed at infrared wave- 
lengths from Earth, consists of three peaks, the first two 
faint and brief and the third bright and long-lasting 
[Graham et ai, 1995; Nicholson et ai f 1995a]. The 
three peaks were seen for all the well-observed events 
(P. Nicholson personal communication, 1995; e.g., see 
Waianabe et ai, [1995], for an excellent recording of 
the bright K event at 2.35 /im), but were first reported 
for the well-observed and well-favored R event. Figure 
1 is a cartoon that illustrates the viewing geometry of 
a typical SL9 event. 

The first peak, or first precursor, preceeded by some 
10 seconds the combination entry flash and fireball seen 
directly by the Galileo spacecraft [Graham et ai, 1995; 
Nicholson et ai, 1995a]. As observed from Earth, the 
first precursor was some 4 orders of magnitude fainter 
than was the entry flash as seen by Galileo (P. Nichol- 
son personal communication, 1995). At first the fireball 
was hot enough that it emitted mostly visible light. But 
as the fireball expanded it cooled. Expansion was effec- 
tively adiabatic [Zahnle and Mac Low, 1994a; Chevalier 
and Sarazin , 1994; Carlson et ai, 1995a]. As the tem- 
perature dropped the emitted radiation reddened and 
the plume faded from visible light over some 30-40 s 
[Chapman et ai, 1995; Martin et ai, 1995]. 
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The second peak corresponds to the fireball itself ris- 
ing high enough to be directly visible from Earth [Gra- 
ham et ai, 1995]. That this would happen was pre- 
dicted by Boslough et ai, [1994] and Ahrens et ai, 
[1994]. By this time (~ 60 s) the fireball had cooled 
to ~500 — 700 K (Galileo near infrared mapping spec- 
trometer (NIMS), [Carlson et ai, 1995a]). Accordingly, 
the second peak was more prominent at longer wave- 
lengths, as is immediately apparent when the Palomar 
light curves (3.2 and 4.5 /im, [Nicholson et ai, 1995a]) 
and the Keck light curve (2.3 /im, [Graham et ai, 1995]) 
for the R impact are compared. The onset of the sec- 
ond peak was abrupt [Graham et ai, 1995]. This can 
be explained by the fireball having a sharply defined 
outer edge at early times, which it has in some numer- 
ical simulations [e.g., Zahnle and Mac Low , 1994a], As 
the fireball continued to expand and cool, it continued 
to fade. 

The fireball expands so enormously that its tempera- 
ture drops below 50 K. At some point along the cooling 
trend it becomes inappropriate to call it a “fireball” . 
Rather, the former fireball becomes the vanguard of the 
ejecta plume. Almost all of its once enormous thermal 
energy is converted to kinetic energy of expansion. As 
it expands first silicates condense, later carbonaceous 
matter, then water (if present), ammonia, and so on. 
The condensates made the sunlit parts of the plume 
visible from Earth, to be imaged by the Hubble space 
telescope (HST) [Hammel et ai, 1995]. 

The third peak corresponds to the plume falling back 
on the atmosphere and the impact site rotating into 
view [Zahnle and Mac Low , 1994b; Nicholson et ai, 
1995a; Graham et ai, 1995]. Interestingly, the timing 
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Figure 1. Cartoon illustrating observing geometries for a typical SL9 impact. 


of the IR light curve recorded by Galileo for the R and 
G events is similar to that observed on Earth [Carlson 
et al. , 1995b]. This implies that something other than 
Jupiter’s rotation delayed the onset of the third peak to 
6 min after impact. In contrast to the first two peaks, 
which account for only a small fraction of the impact 
energy because the fireball has a small surface area, 
thermal radiation produced by the plume’s reentry into 
the atmosphere can account for a large fraction of the 
impact energy, because the plume falls over an enor- 
mous area. Most of this paper addresses the thermal 
radiation produced when the plume falls down. 

Still later bumps in the light curve, the so-called 
bounces, are probably just that, but will not be ad- 
dressed further in this paper. 

Precursors 

Because the impacts occurred on the back side of Jup- 
iter the visibility of the entry flash or fireball on Earth 
needs explanation. One possibility is that terrestrial 


observers saw the fireball itself dimly reflected off high- 
altitude dust, possibly of cometary origin [Hammel et 
al , 1995]. This no longer seems likely, now that it is 
well established that the first precursor preceeded the 
Galileo events by some 10 s. It is likelier that we were 
seeing direct emission from the high-altitude meteor 
trail [Graham et al, 1995; Hammel et al. , 1995]. Iden- 
tification with the high altitude meteor trail is strongly 
favored by the relative timing of the terrestrial and 
Galileo events [Nicholson et al t 1995a]. The energy 
is reasonable. Zahnle and Mac Low's [1994a] nominal 
1 km diameter impactor releases about 0.01% of its en- 
ergy above 250 km above the cloud tops, high enough to 
be directly visible from Earth [Nicholson et al, 1995a]. 
The emission at infrared wavelengths implies a rela- 
tively low temperature. The low temperature argues in 
favor of the meteor trail, effectively a line charge which, 
like the fireball, also expands and cools [Boslough et al, 
1994; Crawford et al, 1994]. 

The fireball rises above the limb about a minute after 
impact. Figure 2 shows fireball light curves calculated 
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Figure 2. Fireball light curves, calculated numeri- 
cally after Zahnit and Mac Low [1994a]. Luminosities 
are extremely sensitive to details in the calculations, 
and should not be taken too literally. The solid curves 
loosely approximate Galileo’s unobstructed view. The 
dotted curves represent the view from Earth. They are 
labeled by the angle past the limb where the impact 
occurred. The scale labeled “Janskys” is an order of 
magnitude estimate only. 

for 10 27 and 10 28 erg impacts. The solid curves show un- 
obstructed light. They loosely approximate what was 
seen by Galileo. The dotted curves show the second 
precursor as viewed from Earth. The labels on the 
dotted curves give the angle 0 O of the impact site in 
degrees beyond the limb. The R impact occurred at 
about Oq =5°. 

These calculations are analogous to the numerical 
light curves predicted by Zahnle and Mac Low [1994a]; 
our method is described there. Initial conditions for 
the models shown in Figure 2 are are described in the 
appendix. The light curves shown here differ from W1 
and W2 shown on Figure 12 of Zahnle and Mac Low 
[1994a] in that they are viewed from the side, they de- 
scribe smaller events, and the wake (a prescribed initial 
condition) is hotter. The double peak shown in Figure 
12 of Zahnit and Mac Low [1994a] represents something 
of a conceptual misunderstanding on our part: our cal- 
culations stitched together two separate models, one for 
the entry and another for the explosion, and we did not 
do this well. 

Figure 2 is meant more as an illustrative cartoon than 
as a serious quantitative model: these calculations were 
predictions, made before the impact, as our answer to 
the predictions of fireball visiblity made by Ahrens el 
ai [1994] and Boslough et al. [1994] We have made 
no effort to account for the wavelength of the observa- 
tion, so that the scale labelled “Janskys” is an order 
of magnitude estimate only. Because radiative cooling 
is only a tiny fraction of the energy budget of the fire- 
ball, calculated luminosities are extremely sensitive to 
details in the numerical model (especially temperature) 
and to modest variations in impact parameters. An 
important limitation is that the only opacity source we 
considered was the H“ ion. The calculated fireball fades 
too quickly because the H _ disappears exponentially 
with cooling temperatures. In the real plume at tem- 
peratures below about 1500 K the main opacity source 


would have been molecules or dust, either silicate or 
carbonaceous. Because it is more opaque, the real fire- 
ball would fade more slowly than the calculated fireball 
[Graham et ai, 1995]. 

Nicholson et ai [1995a] point out that the decay 
times for the second peak are counterintuitive: the 2.3 
micron curve decays much more slowly than the 4.5 
micron curve, although the latter was at first much 
brighter, and further cooling would suggest that fad- 
ing at 4.5 /im would be slower. The resolution may 
be due to different opacity sources (either molecular or 
aerosol), as 2.3 /im is deep in a methane-carbon monox- 
ide absorption feature, and small aerosols emit better 
at shorter wavelengths. 

The Toy Plume 

As a first model for the infrared light curve, we have 
devised a simple two-dimensional (2-D) axisymmetric 
toy model to describe the decline and fall of the plume. 
The toy plume assumes that the velocity field of the 
ejecta plume is ballistic and that the mass-velocity dis- 
tribution in the plume follows a power law that we have 
calibrated to numerical models. Two additional simpli- 
fying assumptions, that the plume is axisymmetric with 
an opening angle 0' , can easily be changed to treat an 
inclined, asymmetric plume. The ballistic approxima- 
tion is valid to a factor of order c 2 <C v 2 , where c $ is 
the sound speed and v a typical velocity in the plume, 
and is in excellent agreement with the results of detailed 
numerical models. 

We divide the ejecta plume by position (r, 0, <£) into 
a vast number of mass elements. Each element is 
launched on its own unique ballistic trajectory. We tile 
Jupiter’s “surface” by distance and azimuth to produce 
a kind of dart board centered on the impact site. We 
then count up where and when the mass elements reen- 
ter the atmosphere (i.e., where and when they hit the 
dart board). The effective radiating temperature of the 
reentry shock is determined by balancing the energy 
supplied by infalling ejecta against thermal radiation by 
opacity sources (dust, soot, molecules, darts, etc.) em- 
bedded in — or generated by — the reentry shock. The 
approach is closely analogous to that used by Melosh et 
ai [1990] and Zahnle [1990] to model thermal radiation 
following the K/T impact. We account for the viewing 
geometry of each surface element when integrating the 
radiated flux as seen on Earth. 

Velocity 

The mass-velocity distribution of ejecta from hyper- 
velocity impacts generally obeys a power law. It is con- 
venient to write the power law in cumulative form, 

M(> u) oc v ~ Q , (1) 

in which the notation M(> v) refers to the cumula- 
tive mass ejected at velocities greater than v. Exam- 
ples abound. In a study of ejecta scaling, Housen et 
ai [1983] cite impact experiments in sand (a = 1.23), 
water (a = 1.66), and basalt (o = 2). The experimen- 
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tally based crater-scaling rules suggested by Schmidt 
and Housen [1987] are equivalent to flow fields in which 
a = 1.33 (nonporous target) or a = 1.8 (porous tar- 
get). The power a is constrained by the conservation 
laws of momentum and energy to be between 1 and 2 
[Zel’dovich and Raizer , 1967, p. 834; Housen ei ai } 
1983]. 

Perhaps more directly relevant to airbursts on Jup- 
iter is an approximate similarity solution for the normal 
impact of a gas “plate” on a gas half-space [ZeVdovich 
and Raizer, 1967, pp. 827-833, in particular, equations 
8.38-8.42 on p. 828]. ZePdovich and Raizer define the 
cumulative mass above an altitude z as m, and the cu- 
mulative mass above the shock front as M . The latter 
is a monotonically increasing function of time. The up- 
ward velocity u increases monotonically with z (in the 
asymptotic limit, u increases linearly with z). An ana- 
lytic solution is obtained for the special case of an ideal 
gas with 7 = 1.4. The cumulative mass m with velocity 
greater than or equal to u is given by 

woe M~ 2/3 ((m/M)~ 2/3 - 3). (2) 

In the limit that m « M, we obtain u oc m~ 2 / 3 . The 
asymptotic mass- velocity distribution of the ejected gas 
therefore takes the form M(> v) oc v~ 1,5 . Zel’dovich 
and Raizer argue that the power a depends on the 
geometry of the blow-out as well as on 7, and so we 
should not necessarily expect a — 1.5 to apply to three- 
dimensional explosions on Jupiter, nor should we expect 
the same value of a to apply to impacts of all magni- 
tudes. 

Figure 3 shows cumulative velocity distributions at 
50 and 100 s calculated from two hydrodynamic simula- 
tions with 7 = 1.4. (see appendix for numerical details). 
The simulations are similar to those published elsewhere 
by Zahnle and Mac Low [1994a]; cumulative velocity 
distributions are presented in their Figure 8). Figure 3 
was constructed by adding up all the mass on the grid 
with velocity greater than v . We ignore the direction 
of motion, although most motion is outwards, as one 



Figure 3. Cumulative mass-velocity distribution for 
two numerical models at times of 50 s (solid curve) and 
100 s (dotted curve) after impact. The dashed line is a 
power law with exponent a = 1,55. 


expects. Immediately apparent from the figure is that 
the cumulative mass distribution does indeed approxi- 
mate a power law for velocities below about 10 km/s. 
The power, a « 1.55, is close to a for the idealized one- 
dimensional impact. This is mostly coincidence, since 
a different numerical model (not shown), featuring a 
relatively shallow explosion, gives ot = 1.35. The distri- 
bution truncates at a maximum velocity, call it r> max , 
above which a much steeper distribution, a ~ 6, ob- 
tains. Two specific cases are shown, one 4 times hotter 
than the other (constructed by setting the mass density 
in the wake 4 times lower). The hotter simulation gives 
higher maximum ejecta velocities. The cumulative ki- 
netic energy in the ejecta plumes is quite large, about 
35-40% of the total impact energy in both numerical 
models. 

The explosion initially produces a small volume of 
very hot gas. This volume expands almost adiabati- 
cally, converting most of its initial thermal energy into 
kinetic energy. Because the firebalFs radiating surface is 
relatively small, radiative losses are small. As the plume 
cools, pressure forces become unimportant and the ve- 
locity field becomes inertial and expansion homologous 
(shape-preserving), characterized by velocities increas- 
ing linearly with distance from the origin, such that 


^ ^max 

R Rmax 


( 3 ) 


where R is the radial distance from the origin, and R m * x 
is the distance to the outer edge of the plume (defined 
by ^maxO- As these trajectories begin to bend under 
the influence of gravity a nearly purely ballistic velocity 
field results. The radial and vertical components of the 
ballistic velocity field are accurately given (rotation is 
ignored) by 


- = 7 (4) 

z 1 /r x 

= J~2 9L (5) 

We have verified this by directly superposing this simple 
velocity field on the velocity field calculated by detailed 
2-D gas dynamical simulations (the 3-D velocity field 
would be more complicated). 


Density 

For the toy plume we assume that the power law 
mass- velocity distribution applies after most of the ther- 
mal energy of the plume has been converted to kinetic 
energy. For simplicity, we also assume that the plume is 
axisymmetric and that it is initially ejected from a cone 
with opening angle 9*\ i.e., ejection is isotropic for all 
angles 0 < 0 * , where the opening angle O' is measured 
from the zenith. Our 2-D numerical models for candi- 
date SL9 events typically explode with ff « 45° [Zahnle 
and Mac Low, 1994a], although this is by no means an 
exact result. In general, the opening angle depends on 
the relative energy of the explosion. A very energetic 
explosion occurring very high in the atmosphere will 
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have a wider opening angle than a low-energy explo- 
sion occurring deep in the atmosphere; the important 
dimensionless discriminant is Ei/p 0 H 3 > where Ei is the 
explosion energy, p 0 the ambient pressure, and H the 
local scale height. Also, the numerical models appear to 
show that the opening angle is larger for high velocity 
material than for low- velocity material. 

We can use the total mass and energy of the plume to 
obtain its density distribution as a function of ejection 
velocity (or, equivalently, as a function of position in 
the plume). Using (3) to relate R and v, the cumulative 
mass M(>v) in the plume is 

M{> u) = 2tt (1 — cos $') ( — — f p(v)v 2 dv. 

\ ^max / J v 

( 6 ) 

It follows from (1) that p(v) oc v 3 The cumulative 
kinetic energy of the plume is 

E(>v) — 7T (1 — cosfl') ( -- - -- ) / p(v)v 4 dv. (7) 

\ ^max J Jv 


Both energy and mass are conserved. If we assume that 
the plume is invested with a fraction r? of the total im- 
pact energy | m,- vf , we obtain the density distribution 


(2 — ot)r)rr\i (t^/umax) 2 /tWc\ 3+a 

2,(1 -co JFP&Z (tv 


( 8 ) 


Equation ( 8 ) is valid for v < i> ma x' In relating ( 8 ) to ( 6 ) 
we implicitly assume that M(> f max ) — 0- Note that, 
because J? max = Vmax^ the density of the plume decays 
as t~ 3 y as expected for a homologous expansion. Our 
numerical results indicate that 0.3 < p < 0.45; i.e., a 
substantial fraction of the impact energy is steered into 
the plume. This value may prove controversial, so we 
will treat 77 as a parameter that others may set to 0.01 
if they are so inclined. 

The normalized form of ( 1 ) is 



Equation (9) has been evaluated for a = 1.55 and 
Umax = 12 km/s, the latter consistent with a plume 
3000 km high. 

When the plume falls back onto the atmosphere, it 
liberates the vertical component of its kinetic energy 
in a reentry shock. There are actually two shocks, one 
driven into the Jovian atmosphere and the other in the 
falling plume. Both the plume and the atmosphere it 
strikes can get quite hot. Peak shock temperatures are 
highest in the more distantly thrown ejecta. The hot 
gas and embedded particulates cool radiatively. 

The horizontal component of the velocity is preserved 
across the shock. Therefore the shocked plume con- 
tinues to expand radially for a considerable time after 


it reenters the atmosphere. The plume expands hori- 
zontally until it has swept up its own mass in Jovian 
air; thus it travels about as far after it reenters the 
atmosphere as it does while in space. The characeris- 
tic crescent-shaped footprint of the plume was observed 
to have rotated through a larger angle than can be ac- 
counted for by Coriolis force acting only while the plume 
was in orbit [Hammtl et al., 1995]. The additional rota- 
tion of the footprint was caused by the plume’s contin- 
uing radial expansion for some 20-30 min after reentry. 

We model the decline and fall of the plume by subdi- 
viding the early plume three dimensionally (r, 9 , <j>) into 
a vast number of mass elements. Each mass element 
is launched on its own ballistic trajectory. We then 
count up where and when the elements reenter the at- 
mosphere. This approach is computationally intensive, 
and hence well-suited to the expensive and astonish- 
ingly powerful machines that would otherwise sit idle 
on our desktops. The mass of a finite element of the 
plume with velocity altitude 9 , and azimuth <f> is 

dm = &j) sin 9 69 p(v)v 2 Sv f (10) 

\ v max J 

The mass element reenters the atmosphere a horizontal 
distance 

r(v, 9) = v 2 sin 29/g (11) 

from the origin at time 

t(v } 9) = 2v cos 9/g } (12) 

as determined by conventional ballistics in the flat Jup- 
iter approximation. 

Pressure 

The pressure level of the reentry shock, if strong, is 



This can be evaluated directly to give p f as a function 
of position and time, but this is not as useful or as 
accurate a measure of “the” shock pressure as one might 
hope. To first approximation, the greatest mass flux 
of material to arrive at any point is the first material 
to get there. This is rigorously true for plumes with 
opening angles 9* < 45°, for which the first material 
to reach a given point is ejected with the lowest total 
velocity of any material to reach that point, and hence 
it is the densest gas to reach that point. Thereafter the 
mass flux declines monotonically and the shock rises 
to progressively lower pressures. Meanwhile the heated 
atmosphere both cools radiatively and also rises to seek 
a new scale height. The accumulated mass fallen at a 
given place (multiplied by g) provides a measure of the 
pressure level of the contact discontinuity between the 
plume and the Jovian atmosphere, and we suggest that 
this is probably a good measure of the reentry shock. 

An illustrative analytic approximation to the pressure 
at the reentry shock can be obtained by replacing v 2 in 
(9) by rg, which gives 
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p(r) = 


dM/dr 
2 tt r 9 
(2 — ot)r} m l 
47T 




7 l-a/2 


r 2 + a/2 


(14) 


When evaluated for a — 1.55 and t> max — 12 km/s, this 
becomes 

. f m i \ /1000 km \ 2 8 

P(r) « 90 t, j ( J pbars. (15) 

Shock pressure is directly proportional to the mass of 
the impactor. Pressures at large distances are quite low. 
At r = 2000 km, a typical shock pressure produced by 
a 10 14 g impact with 77 = 0.4 would be 5 /ibars. 

In Figure 4 we compare the pressure level given by 
(15) (faded solid curve) to results obtained directly from 
our fiducial numerical model (diamonds). This is the 
same numerical model by which we calibrated (9), so 
good agreement might be expected; nevertheless, the 
agreement is gratifyingly good. The diamonds indicate 
where tracer particles swept up in the plume have got- 
ten to after 15 min. The lowermost diamonds trace 
the location of the boundary between Jupiter and the 
ejecta. 


Temperature 


The ordinary Rankine-Hugoniot relations for a strong 
shock require that the infalling material reach a temper- 
ature of 



where m is the mean molecular mass. The origin of 
this equation can also be seen as a balance of the initial 
kinetic energy with the thermal energy c v T of the 
halted gas (for a weak shock, temperatures are higher 
than given by (16)). 

The temperature in (16) is the temperature of the 
gas immediately after it passes through the shock. It 
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Figure 4. The altitude and pressure of the reentry 
shock for a 10 27 erg SL9 impact, 15 min after impact. 
The solid curve shows the location of the reentry shock 
as approximated by (15). The diamonds show tracer 
particles that were swept up by the plume. The lowest 
diamonds trace the boundary between Jupiter and the 
fallen plume. 


can be high. Temperatures would be especially high 
in material that originated from the comet because the 
mean molecular weight of the vaporized comet is rela- 
tively large, probably of the order of 10 < fi < 20. The 
detection of hot (>2000 K) carbon monoxide emission 
at 2.3 fim (R. Knacke, personal communication, 1994); 
[Meadows ct ai , 1994] offers direct evidence that for- 
mer cometary material became very hot on reentering 
the atmosphere. 

The effective radiating temperature, between 500 and 
1000 K [Nicholson et ai , 1995b], was much lower than 
implied by (16). Evidently radiative cooling by the 
shocked material was efficient. HST images clearly re- 
vealed that the ejecta blankets were dusty [Hammel et 
ai, 1995]. The dust is an obvious opacity source for 
radiative cooling. 

Perhaps the simplest expression for the radiating 
temperature of a dusty layer heated by falling ejecta 
is to assume an instantaneous energy balance between 
the kinetic energy of falling matter and radiative cool- 
ing, buffered by the heat capacity of the falling matter 
[Zahnle, 1990; Melosh et ai, 1990]: 

e=— — m + 2oT 4 (1 - e~ 2r ) . (17) 

7 — 1 m 

The factor 2 crT 4 assumes a grey radiator that emits 
both up and down. The factor 1 — e~ 2r allows for the 
optical depth r of the radiating layer; the factor of 2 
gives the right limit as r — ► 0. We solve (17) for the 
temperature T by Newton’s method. 


Opacity 


Radiative cooling depends on the presence of an opac- 
ity source, presumably particulates either in the plume 
or generated when the plume strikes the atmosphere. 
The latter requires a sufficiently strong shock. Zahnle 
et ai [1995] show that reentry velocities must exceed 
4-5 km/s for shock heating to be strong enough to gen- 
erate carbonaceous dust from Jovian air. This can ex- 
plain the relative transparency of the inner parts of the 
ejecta blanket as observed with the HST [Hammel et ai, 
1995], and it is probably required to explain the late 
onset of infrared radiation as observed by the Galileo 
spacecraft [Carlson et ai, 1995a, b]. The observed par- 
ticulates are apparently rather small. Nicholson et ai 
[1995b] suggest an effective optical depth of about 0.022 
at A 0 = 4.5 /im, with a A~* dependence. Molecular 
opacity might also contribute. Here we will use a sim- 
ple grey (wavelength-independent) opacity. 

We assume that r is proportional to the total inte- 
grated mass of material falling at radius r, 


r(r, 



m dt ; 


(18) 


i.e., the opacity sources are embedded uniformly in the 
shocked air, and all the dust that falls at any point 
is available to contribute to radiative cooling. This ap- 
proximation ignores the vertical structure of the reentry 
shock (late arriving material shocks at higher altitudes) 
and radial gradients in the plume’s composition (the 
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more distantly thrown material may be comet-rich, and 
hence carry more grains). 

Figure 5 shows effective radiating temperatures at 5, 
10, and 15 min, as a function of distance from the im- 
pact site, calculated with a constant k = 100 cm 2 /g 
(solid curves). This is the plume used to generate the 
light curves shown in Figures 6 and 7. The temper- 
ature field forms warm rings (seen nicely at 10 min). 
The inner regions are cooler because the optical depth is 
relatively high (from early-arriving ejecta) and because 
temperature is calculated by the instantaneous energy 
balance between e and radiative cooling. If either of 
these assumptions is relaxed the expected interior tem- 
perature warms. To illustrate this, the dotted tempera- 
tures on Figure 5 were generated from the same model, 
modified to limit optical depth to material arriving in 
the previous 90 s: 

r(r, t) = k f mdt. (19) 

Jt~ 90 

The temperatures that result are much more nearly con- 
stant, both in radius and time. 

Geometry 


The impact occurs at an angle 0 o over the horizon. 
A given element of the plume reenters the atmosphere 
a time t later, at a distance r from the impact site, 
and at an azimuth <f> (measured from north). During 
this time the impact site rotates into view at a rate of 
0.01° s -1 . The projected surface area of the plume’s 
footprint as viewed from Earth is proportional to the 
cosine of the angle ip between the surface normal and 
the line to Earth 


cos \p « sin 


r sin <p 
Rj cos /, 


0 O + 1.77 x 1(T 4 <) > 0, 

( 20 ) 


where « 45° denotes the latitude of the impact, and 
Rj refers to Jupiter’s radius. Where cosxp < 0 the 
surface element is on the far side of the planet. The 



Figure 5. Effective radiating temperatures at 5, 10, 
and 15 min produced by the toy plume with constant 
opacity k = 100 cm 2 /g. The solid curves use (18), to 
which all the particles contribute to optical depth. This 
is our nominal standard model. The dotted curves use 
(19), to which only the most recently arriving particles 
contribute to the effective optical depth. 



Figure 6. Light curves at selected infrared wavelengths 
generated by the toy plume model, as seen from Earth. 
The Keck light curve [Graham et ai, 1995] at 2.3 /im is 
shown for comparison (solid curve). The dashed curves 
(standard model) uses (18) for r(r, t)\ the impactor is 
500 m diameter and 0.5 g/cm 3 . The dotted curves use 
(19) for r(r, t) and a 450 m impactor. 


flux at Earth is 

F v - cos ipda (1 - exp (- r v sec ip )) , (21) 

where B v is the Planck function, A the Earth-Jupiter 
distance, and da a surface element on Jupiter. 

We have neglected azimuthal asymmetry in (21); this 
will be important to a complete model. We have also 
neglected the altitude of the radiating layer. Altitude 
can strongly affect visibility. When we attempt to ac- 
count for this, by using the numerical equivalent of (14) 
and (15) for pressure, visibility from Earth generally 
improves. The result is a brighter, earlier flash. For the 
present we choose to omit this parameter. In a future 
study, one that includes a nonaxisymetric 3-D plume, 
we will consider the role of altitude more fully. 

Light 

Figure 6 shows calculated light curves at 2.3, 3.5, and 
4.5 /im, as seen from Earth, generated by the toy plume 
models (dashed, dotted lines). The dashed curves use 
(18) to prescribe optical depth. The Keck light curve for 
fragment R [Graham et ai, 1995] at 2.3 /im is shown 
for comparison. The match at 3.5 /im and 4.5 /im is 
comparably good (or bad, depending on one’s taste in 
these things; see Nicholson et ai [1995a]. This model 
uses a 500 m diameter impactor of density 0.5 g/cm 3 , 
7] = 0.5, u m ax = 12 km/s, and an opening angle 6 ( = 
45°. The total impact energy is 6 x 10 26 ergs. The value 
of u ma x> 12 km/s, implies a plume that reaches 3000 km 
high. The value of f ma x is set by the duration of peak 
light, and hence is not very negotiable. 

We obtain a good match to the relative fluxes at 2.3, 
3.5, and 4.5 /im using a constant opacity of k — 100 
cm 2 /g, which is fairly high. If we assume particles of 
radius 0.25 /im, density of 1 g/cm 3 , and assume that 
opacity is equal to the sum of their geometric cross sec- 
tions, the particle density in the ejecta plume would 
be 0.004 g/g. If instead we assume that their effective 
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Figure 7. Light curves at selected infrared wavelengths 
generated as seen from directly above the impact site. 
This view more closely approximates Galileo’s. The 
dashed curves are from our standard model. The on- 
set is too early. The dotted curves restrict dust opacity 
to air experiencing shock velocities exceeding 4.5 km/s. 
This approximates the minimum shock heating required 
to generate carbonaceous particulates from Jovian air 
[Zahnle et al. } 1995]. The “spikes” in the light curves 
are probably artificial, caused by relatively high tem- 
peratures when optical depth is very low. 


cross sections are proportional to A -1 , the required par- 
ticle density in the plume rises to 0.04 — 0.08 g/g. Such 
a high density would require that a fair fraction of the 
plume is of cometary origin. According to Pollack et al. 
[1994], the Eosseland mean opacity of a solar composi- 
tion mix of gas and dust, in which 60% of the carbon is 
in organic grains, is about k — 1 cm 2 /g; for shocked Jo- 
vian air the corresponding opacity would be k « 20 — 30 
cm 2 /g. Again, to raise the opacity to k « 100 cm 2 /g 
requires a significant admixture of comet. 

The calculated light curves using (18) for r (dashed 
lines) decay too abruptly, especially at 2.3 pm. This oc- 
curs in our standard model because (1) the temperature 
drops rapidly once most of the mass of the plume has 
reentered the atmosphere, caused by the optical depth 
(18) remaining high while e declines, and (2) the plume 
is artificially truncated at v < i> max , which effectively 
terminates the event at t = 2t> max /</. There are several 
effects that could act to smooth the light curve. The 
continued radial spreading of the ejecta pattern after 
reentry, alluded to above, is one effect. A second is fi- 
nite heat capacity. We have assumed an instantaneous 
energy balance. Allowing for history would moderate 
the temperature decline. A third, also discussed above, 
is that the optical depth of the earliest arriving ejecta 
may not be relevant to the late arriving ejecta. On 
the other hand, it is interesting to note that the gen- 
eral form of the computed light curve at 2.3 /im is very 
similar to that of the bright K event [Watanabe et al. , 
1995]. It would appear that the smaller, optically thin- 
ner events produced smoother light curves, while the 
larger, optically thicker events were more asymmetric. 

The dotted curves on Figure 6 show how the light 
curves change when (19) is used to prescribe r. The re- 
sult is that the light curves at 2.3, 3.5, and 4.5 /im are 


smoother and more nearly parallel, as observed [Nichol- 
son et al . , 1995a]. To match the observed fluxes this 
model uses a 450 m rather than a 500 m impactor. 
In this model t> max = H km/s (the best fit is with 
Umax = 10.8km/s). 

Figure 7 shows the same event (450 m object) as seen 
from directly above the impact site. These light curves 
approximate the view from Galileo. Because the radi- 
ating layer is mostly optically thin, the fluxes seen from 
Galileo are not likely to be much smaller. The event 
is brighter when seen centered on Jupiter’s disk, and it 
starts earlier. The predicted early rise (dashed curve) 
was not seen by Galileo for the R and G impacts [Carl- 
son et ai, 1995b]. The delayed onset of radiation from 
the reentry of the plume is not plausibly explained by 
an absence of low-velocity ejecta. A better explana- 
tion, and one that receives direct visible support from 
the spectacular I1ST photos, is that low-velocity ejecta 
were largely transparent. Low-velocity ejecta are prob- 
ably relatively comet-free, hence to become opaque the 
air must generate its own opacity. Zahnle et ai [1995] 
show that if reentry velocities exceed 4.5 km/s the shock 
is strong enough to generate carbonaceous dust, possi- 
bly soot, from Jovian air. The dotted curves show the 
consequences on the Galileo light curve of restricting 
dust formation to air shocked at v z > 4.5 km/s. The 
fast rise at 6 min is a probably artificial consequence 
of turning opacity on abruptly at 4.5 km/s; when the 
opacity is first turned on, the temperature is too high. 

Taken at face value, the toy plume model implies that 
the R event was produced by the impact of a 450 to 500 
m diameter body of mass 2 x 10 13 to 3 x 10 13 g. By ex- 
tension the largest plumes would be created by 1 x 10 14 
to 2 x 10 14 g objects. These sizes are consistent with the 
tidal-breakup models of Asphaug and Benz [1994], and 
Solem [1994]. However, energy estimates made by the 
light of the plume are somewhat sensitive to opacity and 
viewing geometry. The view from Galileo is less sensi- 
tive to geometrical details and hence should provide a 
better observational constraint on the actual energy of 
the plume. 

Appendix: Initial Conditions to 
Numerical Models 


The initial conditions of the explosions assume the 
bolide wake can be represented by a moving line charge 
with an energy deposition rate described by the ana- 
lytic theory of Zahnle and Mac Low [1994a]. There 
we made the assumption that the impactor expands 
quasi-statically due to the pressure difference between 
its front and sides to derive an approximate analytic 
expression for its downward velocity, 


v — Vi exp 


f — nC^piz) 2 H 3 sec 3 8 

K 2 tm pi 


(Al) 


where the impactor has an initial downward velocity 
at an angle to the horizontal of 6 , a density of pi, and 
a coefficient of drag C/j, which has values of 1.0 for a 
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sphere and 1.7 for a cylinder. The atmosphere has a 
density profile p and a scale height H . 

We need to make one further approximation in or- 
der to derive an analytic solution for the energy release 
profile. We approximate the density profile as a simple 
exponential, p = po exp ( — z/H). In the simulation we 
actually use a more realistic atmosphere with variable 
scale height [see Zahnle and Mac Low, 1994a], so we 
take H to be the value near the point of maximum en- 
ergy release. We also add a term to account for high 
altitude radiative ablation to get 


1 dE 
p dz 


2 pCj)H 


+ -Cp nrfv 2 sec 6 , (A2) 


where the initial radius of the impactor is r;. 

The radius of the wake is a parabola in altitude z , 


1/2 1 dE_ z_--z 1 
p dz v,- ’ 

(A3) 

where z 0 is the airburst altitude [Zahnle and Mac Low, 
1994a; equation 5]; r* corresponds to the radius of the 
impactor, and Vi to its velocity. We place the bottom 
of the wake at an altitude z f slightly below z 0 , at 

z' = z 0 - (2 C D p/pi) 1/2 H seed - r,. (A4) 



The wake is assumed to be well-mixed so that energy 
and mass are uniformly distributed across the wake. Its 
velocity is also neglected, since a moving line charge has 
no momentum. The energy density in the wake is then 


e'(z) = e(z) + 


1 dE_ 
nR 2 (z) dz 


(A5) 


and the mass density in the wake is 

, 1 dm 

^ ^ 7 rR 2 (z) dz 

The assumption of a uniform wake is not a very good 
one. In reality most of the mass and energy of the wake 
are concentrated at the surface of the paraboloid, while 
the interior is mostly evacuated. However, the numer- 
ical model requires several grid elements to resolve a 
shock. In practice more realistic initial conditions re- 
quires impractically high resolution in the wake. 

The models used to generate the light curves in Fig- 
ure 2 take z 0 = —80 km (10 28 ergs, r* = 0.5 km) and 
z 0 = —25 km (10 27 ergs, r* = 0.236 km); both use 
7 = 1.2. The models used to generate the cumulative 
mass- velocity curves in Figure 3 take z 0 = —25 km, 
r x = 0.297 km, m,- = 1.1 x 10 14 g, E{ = 2 x 10 27 ergs, 
and 7 = 1.4; the hotter model (larger r max ) was gen- 
erated by reducing the mass density in the wake by a 
factor of 6.8. 
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